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2.2 Lubrication approximation for flow in a thin layer 

An essential first step of any analytical approximtation is the art of scaling, which we shall 
emphasize repeatedly throughout this course. 

Let H be the characteristic depth and L the characteristic length in the direction of the 
flow, and assume a shallow layer, i.e., 

if/L«l (2.2.1) 

Let U be the scale of u, then by continuity, the scale of v must be Uj^ in order not to violate 
mass conservation. Leaving the velocity and pressure scales U, P undermined for the time 
being, we introduce the following scales and normalized variables, denoted by primes, 

TT 

t = Tt', x = Lx', y = Hy', u = Uu' , v = U—v', p = Pp , (2.2.2) 

Lj 

The normalized continuity equation is 

Both terms are equally important, reflecting the holyness of the law of mass concentration. 
The longitudinal momentum equation is normalized to 

Udu' U 2 ( ,du' ,du'\ n P dp' vU (H 2 d 2 u' d 2 v'\ , 

TW + T{ u ^ + v W) =gsm9 ~^ + lPWd^ + dy^) {22A) 

Dividing by we get 

H 2 du UHHf ,du' ,du'\ _ gsm6H 2 PH 2 dp' ( H 2 d 2 u' d 2 u'\ 

^~dl + ~^T{ U dx^ + V dy^)~ vU P LvUlh J+ {jJdx^ + dy^) (2 ' 2 ' 5 ^ 

For a shallow layer (H/L 1) we assume in addition, 

UH 



v 

and 



= 0(1) (2.2.6) 



^ « 1 (2-2.7) 
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Omitting terms of the order H/L and smaller, the above equation can be approximated to 
the leading order by 

gsmdH 2 PH 2 dp' d 2 u' 
~ vU pLuUdx^ + dy l2 ( ' ' ' 

or in dimensional form, 

= 9 sin S -I^ + ! £H (2.2.9) 
p dx dy 1 

All inertia terms are inconsequential; the most important balance is among gravity, the 
pressure gredient and the dominant viscous stress. This balance also implies a pressure 
scale, 

P=^L ,2.2.10) 
From the transverse momentum equation, 



H 
T 



Udv' U 2 ( ,dv' ,dv'\ 

1 U h V 

Tdt' L \ dx' dy'J 



n P dp' HuU d 2 v' d 2 v'\ , nn „ s 
9^9-— 7 JL + + (2.2.11) 



pH dy' L H 2 \ dx' 2 dy 



or 



H 2 (H 2 dv' UHH( ,dv' ,dv r 
~T? \yf~dV + \ U dx' +V dy' 



gsinOH 2 H PH 2 dp' H 2 ( H 2 d 2 v' d 2 v r 



vU Lt&n9 pLvU dy' L 2 \L 2 dx' 2 dy' 2 
Either for finite bed slope or for small slope but 



(2.2.12) 



O^)=tan0«l (2.2.13) 

the left hand side above is negligible with an error of 0(H/L) 3 . In physical variables the 
approximate result is 

= -£cos#-i^ (2.2.14) 
pdy 

Not only the inertia terms are insignificant, the pressure is hydrostatic. This balance also 
implies the pressure scale 

P = pgHcos9 (2.2.15) 
Note that (2.2.10) and (2.2.15) together implies the velocity scale 

H qH cos 9 

U = —- (2.2.16) 

The distinguishing feature of negligible inertia is shared by the slow flow through thin 
gaps of bearings in the theory of lubrication. Hence (2.2.9) and (2.2.14) can be called the 
lubrication approximation. 
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We leave it as an exercise to show by similar normalization, that the dynamic boundary 
conditions on y = h can be approximated to the leading order by 



for the tangential stress, and 



du 
dy 



p = 



for the normal stress. It follows by integrating (2.2.14)) that 

p(x, y, t) = pg cos9[h(x, t) - y\ 
The longitudinal momentum equation can also be readily integrated, 



P9 ( . , „dh\ (y 2 ' 

u = sin — cos 6— hy 

H V dx I V 2 



The total discharge is 



fh 

Q = uh = udy 
Jo 

which can be inserted in (??) to give 



pgh 3 ( . dh' 

sin V — cos U— 

3/x V ox 



dh pg d 
dt 3p dx 




This is a nonlinear diffusion equation governing the evolution of the fluid depth. 
In the special limit of a uniform flow, d/dx = 0. The velocity profile is then 

pgh 2 fy y 2 \ 

u — sin v - — — — 

H \h 2h 2 J 

with h being a pure constant. The coresponding discharge is 

Q = — — sin 
Sp 
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